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Abstract 

We obtain a general characterization of discrete-time all-pass rational matrix functions from state- 
space representations. It can be employed to address model reduction problems in the same vein of the 
theory developed by Glover in the continuous-time. Besides model reduction, this characterization is 
shown to be useful in a variety of contexts such as studying LMI’s and Riccati equations and especially 
in the factorization of all-pass functions. The results are obtained in the most general setting, without 
introducing any ad hoc assumption. 


I. Introduction 

In this paper we provide a eompletely general eharaeterization and parameterization of diserete- 
time all-pass matrix funetions and use this result to deseribe in full generality the geometry of 
the solution set of eertain LMI’s and of the assoeiated Rieeati equations. We also develop a 
faetorization theory and related state-spaee proeedures for the faetorization of all-pass funetions. 
The eharaeterization of diserete-time all-pass matrix funetions presented in the main theorem 
of the next seetion, parallels the eontinuous-time fundamental result of Glover’s flU Theorem 
5.1] in the most general setting, without introdueing any ad hoc assumption. A detailed rigorous 
proof of this result seems to be presented here for the first time after past unfruitful attempts in 
the literature. Commonly used faeilitating assumptions in diserete-time, sueh as non-singularity 
of various matriees (in partieular of the A matrix) and unmixing are avoided and only diseussed 
as partieular eorollaries of more general statements. 

A diserete-time version of Glover’s model reduetion proeedure seems to be worth as diserete- 
time models are often the rule in applieations. Its derivation however is not a simple transposition 
of the arguments used in eontinuous-time as there are several differenees whieh make the job 
teohnieally mueh harder, see for example the attempts in lUOl . Some of these diffieulties are 
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well-known. The use of the often advocated Cayley transform requires for example invertibility 
assumptions which are not met in some applications (see, e.g., the comment in [fTTl p. 1996]). 
Moreover, it seems to be an accepted point in the systems and control community that, as stated 
in IflSl p. 559], “... it is generally more appealing to give derivations in the coordinates of the 
original [discrete-tme] data; also algorithms may be more reliable if generated for the specific 
model class”. Apparently a discrete version of the continuous-time all-pass dilation of Glover 
under general hypotheses as those made in the present paper has been lacking. So far, to our 
best knowledge, the book literature of the last two or three decades, e.g. [[D or, ifT^ flUl seems to 
be just re-proposing continuous-time model reduction and does not address a discrete-time 
version of Glover’s theory. 

The results of the paper have many possible applications. Applications to Hankel-norm ap¬ 
proximation of rational discrete-time transfer functions may now be pursued by just following 
the route shown in the paper In Chapter 16 of the book ffTTII a slightly less general character¬ 
ization of discrete all-pass functions is used to do Hankel-norm stochastic model approximation. 
Stochastic modeling without stability constraints is another direction which has been touched 
upon in [|6l, further exposed in [fTTl and can be addressed in wider generality by using the 
techniques described in this paper. This is a relatively unappreciated area of stochastic modeling 
which has several applications to smoothing and to non causal estimation. We believe that this 
setting is worth understanding especially because of a very illuminating isomorphism with LQ 
control with an indefinite cost function. In a companion paper we shall apply this isomorphism 
to resolve an old open problem about the existence of negative semidefinite solutions of the 
Riccati equation of LQ control. 

The lay-out of this paper is as follows: 

Section |I^ contains the statement and proof of the main result. The proof is essentially self- 
contained save for a technical Lemma from [[3l which considerably generalizes a result on 
controllability due to Wimmer lfT5]l . 


In section III we introduce two dual linear matrix inequalities with a rank constraint which 
define families of square all-pass functions having a fixed pole structure. We prove a geometric 
characterization of all solutions in terms of A- or A'^- invariant subspaces. When A is non 
singular these matrix inequalities turn into two dual homogeneous algebraic Riccati equations. 
A very exhaustive classification and description of the solutions of those Riccati equations is 
provided. It is well-known, see e.g. IfTTlI that the analysis of algebraic Riccati equations can be 
reduced to that of homogeneous Riccati equations. 
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The study of families of solutions of the constrained LMI’s of Section III unveils the basic 
principles and a direct method to characterize and classify the left- and right all pass factors of an 
arbitrary square all pass rational function. Rational factorization theory was first systematically 
discussed in the early book [|2l| quite heavily relying on the assumption of an invertible D 
matrix. Here we extend the factorization results of Fuhrmann and Hoffmann [|71 derived for 
inner functions, under general hypotheses. When A is non-singular the classification can be 
given directly in terms of solutions of two dual homogeneous algebraic Riccati equations. 

In the concluding section we indicate some possible generalizations to non square matrix func¬ 
tions. 

Notations in the paper are quite standard; we only mention that denotes the MoorePenrose 
pseudoinverse of the matrix X. A technical condition which is often referred to is that of 
unmixing. One says that A G has unmixed spectrum or, briefly, is unmixed if it does not 

have reciprocal pairs of eigenvalues. In particular an unmixed matrix cannot have eigenvalues 
of modulus one. 


H. The main result 

Theorem 2.1: 

1) Let 

Q{z):=C{zI -A)-^B + D (1) 

be a minimal realization of an m x m rational discrete-time all-pass function. Then A is 
non-singular if and only if D is non-singular. 

2) Let Q be a minimal realization of a rational discrete-time all-pass function. Then there 
exist two invertible matrices P = and Q = such that PQ = I and 

APA~^ -P = BB'^ 

< BD^ - APC^ = 0 ( 2 ) 

DD~^ - CPC^ = I 

A~^QA -Q = C^C 

< C^D - A~^QB = 0 (3) 

D^D - B^QB = I 

3) Let Q be a minimal realization of a rational discrete-time all-pass function. If equations 
Q admit a solution P, then such a P is unique. If equations Q admit a solution Q, then 
such a Q is unique. 


r\r> A TT'X' 
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4) Let A G B G C G D G be given (no minimality is now 

assumed). If there exists P = satisfying (|^ then Q{z) given by ([^ is all-pass. 

Similarly, if there exists Q = satisfying (|^ then Q{z) given by ([^ is all-pass. 

Finally, P is a non-singular solution of (|^ if and only if P~^ is a (non-singular) solution 
of©. 

5) Let A G B G be a given reachable pair. Then, P = P^ is such that 

APA~^ - P = BB'^ (4) 

if and only if there exist matrices C G and D G such that Q{z) given by 

Q is a minimal realization of an all-pass function and P is the solution of ([^ for the 
quadruple {A, B,C, D). In this case, P is necessarily non-singular and such that 

I + B^P-^B >0. (5) 

6) Let A G C G be a given observable pair. Then, Q = is such that 

A~^QA -Q = C^C (6) 

if and only if there exist matrices B G and D G such that Q(z) given by 

Q is a minimal realization of an all-pass function and Q is the solution of Q for the 
quadruple {A, B,C, D). In this case, Q is necessarily non-singular and such that 

I + CQ-^C^ > 0. (7) 

7) Let A G B G C G be given. If there exists P = P^ and Q = such 

that 

APA~^ -P = BB'^ 

< A'^QA -Q = C^C (8) 

PQ = I 

V 

then there exists a matrix D G such that Q{z) given by Q is all-pass. 

Proof: 

1) By assumption we have 

Q{z)Q*{z) = L (9) 

Notice that Q(oo) = D so that by taking the limit 2 ; —)■ cx) in (j^, we see that D is non-singular 
if and only if Q*{z) is bounded at infinity or, equivalently, if and only if Q{z) is bounded in 
a neighborhood of the origin. By taking into account that Q is a minimal realization, this is 
equivalent to A being non-singular. 


r\r> A TT'X' 
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2) Let us first assume that D is non-singular. By recalling point[^, we have that A is non-singular 
as well. 

We have the following minimal realizations: 

Q{z)-^ = D-^ - D-^C{zI - T)-^BD-\ T := A- BD-^C. (10) 

and 


Q*{z) = B'^{z-H - A'^y^C'^ + D'^ 

= d 1 - B'^A-^{zI - A-^y^A-^C'^, (11) 

with Dq := — B^A~^C^, so that, by imposing Q{zy^ = Q*{z), we conclude that there 

exists a unique invertible matrix T such that 

T-^A-^T = A - BD-^C{=T) (12a) 

T-^A-^C^ = BD-^ (12b) 

B'^A-^T = D-^C (12c) 

D-^ = - B'^A-^C'^ (12d) 


By inserting (12c) in (12a) and multiplying on the right side by {A ^T) we get 

T-i _ ^ 


(13) 


so that the first of ([^ admits a solution P = T~^. Moreover, by inserting the expression of 
D~^ provided by (12d) in (12b) we get BD^ = (T“^ + BB'^)A~^C'^, which, in view of ( [T^ , 
may be written as BD^ = AT~^C^, so that P = T~^ solves also the second of ([^. Finally, 
by multiplying (12d) on the left side by D and taking into account of (12c), we easily see that 
P = T~^ solves also the third of ([^. Similarly we see that from (12) it follows that T solves the 
three equations (|^. The proof that T is symmetric is a bit lengthy and is deferred to Appendix 


So far we have established our result in the case when D is non-singular. We now show how 
this case may be viewed as a first step for proving the result in the general setting in which D may 
be singular. Consider an arbitrary rational proper all pass function Q{z) and the corresponding 


factorization (79) established in Lemma A.l of the appendix. Let Qo{z) := Co{zI—Aoy^Bo+Do 
be a minimal realization of Qoiz) so that Dq = Qoioo) is non-singular. Then equations ([^ 
with A = Aq, B = Bq, C = Cq and D = Dq have a symmetric solution Pq which is non¬ 


singular. In view of Lemma A.2 we know that Qi{z) := Qoiz)Qi{z) has the reachable realization 


r\r> \ TT'X' 
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,Si: = 


Qi{z) = Cl {zl — Ai) ^ Bi + Di where 

Cl := [Dq 2 I Co])^i •= 

and Di := [Do,i I 0 ]? 7 i. Now it is immediate to eheek by inspeetion that 

Pi-.= 


0 0 
-Bo ,2 Aq 


0 I 
-Bo,i 0 


Ui, 


I 0 
0 Po 


(14) 


solves equations Q with A = Ai, B = Bi, C = Ci and D = Di. We ean iteratively repeat 
this argument for Qi{z), i = 2,3,..., k and eventually find that Q(z) has a reachable realization 
Q{z) = C [zI — ^) ^ B + D and that equations (j^ with A = A, B = B,C = C and D = D, 
have a solution P. Without loss of generality we may assume that A, B, C are in the Kalman 
reachability form 




A 

0 


B 

c = |C 10 ], 

A = 

A 21 

^22 

, P: = 

B 2 

and P is partitioned conformably as 


p 

P12 




P = 






-'12 

P22 




(15) 


(16) 

By writing block-wise equations ([^ with A = A, B = B,C = C,D = D, and P = P we see 
that the (1,1) block P is a symmetric solution of equations (|^ with A = A, B = B,C = C, 
D = D corresponding to the minimal realization 

Q{z) = c(zI-A) ^B + D. (17) 


The original minimal realization of Q{z) is related to ( [Tt] ) by a change of basis so that there exists 
a non-singular matrix T such that A = T~^At, B = T~^B, C = CT. Then it is immediate to 
check that P ;= is a solution of equations (|^ for the original realization ([^. Observe 

that by minimality of the realization {A, B, C) we have that P, solving the Lyapunov equation 
in Q, is non-singular. 

By resorting to a dual argument we establish the existence of a non-singular matrix Q = 
solving Q. 

It remains to show that PQ = J. To this aim, write (|^ in the form 


PXP^ = X 


( 18 ) 
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where 


Clearly, X is non-singular and 



A 

B 


P 0 

P : = 

C 

D 

X := 

0 -I 


(19) 




0 


( 20 ) 


p-i 

0 -I 

Thus FXF' X~^ = L Therefore, F is non-singular as well and we have F~^X~^ = X~^F~^ 
or 

F^X-^F = X-\ (21) 

The expression (^) of X~^ implies that P~^ is a solution of equations ([s]). As proven below 
these equations however admit a unique solution so that P~^ = Q or, equivalently, PQ = L 
3) Assume that Pi and P 2 are solutions of @ and let A ;= Pi — P 2 . We need to show that 
A = 0. It is immediate to eheek that A satisfies the equations 

^ AAA^- A = 0 

AAC^ = 0 (22) 

CAC^ =0 

From the seeond and the third of these equations we see that im (AC^) is eontained in the non¬ 
observability subspaee of (A, C) and, sinee (A, C) is assumed to be observable, this means that 
AC^ = 0. This implies that CA = 0 and, in turn, C'A^A(A^)^ = 0 for all fc = 0,1,..., n —1, so 
that CA^A(A^)" = 0 for all /c = 0,1,.. ., n — 1. This means that im (A(A^)") is eontained in the 
non-observability subspaee of (A,C), so that, as before, A(A^)” = 0. Now, by multiplying the 


first of (22) on the right side by (A^)”“^, we get A(A^)’^“^ = 0 and, iteratively, A(A^)’^“^ = 0, 
and so on, up to A = 0. The proof for equations ([^ is dual and is therefore skipped. 

4) Assume that equations (|^ admit a solution P = P^. Let us eompute the produet 

$ := Q{z)Q^{z-^) 

= [C{zl - A)-^B + D][B^{z-H - A^)-^^^ + D~^]. 

The first of equations (|^ ean be rewritten as 

BB~^ = {zl - A)P{z-^I - A'^) - zP{z-^I - A~^) 

-z-\zI-A)P, 

so that 

C{zl - A)-^BB^{z-^P 


A 


Ti-ir^T 


C = 


r\r> \ TT'X' 












CPC^ - zC{zI - A)-^PC^ - z-^CP{z-^I - A^)-^C^. 


Moreover, from 


z{zl - A)-^ =1 + A{zl - A)-^ = I+{zI- A)-^A 

it follows that 

C{zl - A)-^BB^{z-H - A'^y^C^ = 

-CPC^ - C{zl - Ay^APC^ - CPA^{z-y - A'^y^C^. 

In eonelusion, we have 

Q{z)Q^{z-^) = DD^-CPC^ 

+cyi - Ay\BD^ - APC^) 

+{DB~^ - CPA'^){z-y - A^y^C'^. 

By taking into account the second and the third of equations ([^, we now get Q{z)Q^{z~^) = I. 

Assume now that equations (|^ admit a solution Q = . By computing the product {z~^)Q{z) 

and using the dual of the previous argument, we get {z~^)Q{z) = I. 

The fact that P is a non-singular solution of ([^ if and only if P~^ is a non-singular solution 
of Q can be shown by defining F and X as in ( [T9l ) and using the same argument that led to 

5) One direction is an immediate consequence of point [^. For the converse, since (A, B) is, by 
assumption, reachable, the solution P of Q is invertible [|3l Lemma 3.1]. Let (n+,n — n+,0) 

-Q 0 

be the inertia of P which is equal to the inertia of Q := P~^. Let E : = 


The inertia of E is given by the inertia of 

__ 0 Im 

of the corresponding Schur complement S which is given by 


-Q 0 


0 Im B'^ 

A B -Q-i 

, i.e. (m + n — n+, 0), plus the inertia 


1 

o 

1 

-1 

A^' 

0 Im 


pT 


^ := -Q-^-[A B] 

= —P + AP AJ — BB^ = Ofixn- 

In conclusion, the inertia of E is (m + n — n+,77,+ ,?7,). On the other hand, the inertia of E is also 
given by the inertia of —Q~^ = —P, i.e. (n — n+,n+,0), plus the inertia of the corresponding 
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Schur complement W whieh is given by 


W := 


-Q 0 




i-Q-^)-^[A B] 


A^QA - Q A~^QB 
B^QA B^QB + I 


(23) 


Henee the inertia of W is given by the inertia of E, i.e. (m + n —n+,n+,n) minus the inertia of 
—Q~^ = —P, i.e. (n — ?7,+ , n+, 0), whieh amounts to (m,0,n). Thus, W is positive semidefinite 
and has rank equal to m. Therefore, there exists a full row-rank matrix [C \ D] G ]^mx(n+m) 
sueh that W = [C \ DY[C \ D]. This means that for the given A and B and for the C and D 
obtained by previous developments, there exists a. Q = P~^ solving ([s]). Then, in view of point 
1^, the eorresponding Q{z) given by Q is all-pass. 

We now prove ([^. Indeed, we have already proved that P~^ solves ([s]) and from the third of 
these equations Q follows immediately. 

It remains to show that {A, C) is an observable pair. To address this issue we exploit ([^ 
whose validity we have already proven. Assume now by eontradietion that the pair (A, C) is 
not observable and let V he a full eolumn-rank matrix whose eolumns (at least one by the 


eontradietion assumption) form a basis for the unobservable subspaee := ker 


C 

CA 


so 


that 


CA^-^ 


CV = 0 


(24) 


and sinee is A-invariant, there exists a matrix K sueh that 


AV = VK. (25) 

By multiplying the first of ([^ on the right side by V we get /C' QAV = QV. We now multiply 
the first of (|^ on the right side by AV and on the left side by A^: We get {A'^YQA'^V = 
A~^QAV = QV. We ean iterate this argument and multiply the first of (|^ on the right side by 
A^V and on the left side by (A^)*^, fc = 2, 3,..., getting 

{A~^yQA^V = QV, 1 = 1 , 2 ,.... (26) 

We now show that 

U := QA^V Y 0, (27) 


r\r> A TT'X' 
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where n is the dimension of A. In faet, from (26) we get {A^)'^U = {A~^)'^QA'^V = QV and 
sinee Q is non-singular and V has full eolumn-rank this yields ( |27| ). We now eonsider the seeond 
of equations ([^. From this equation, we get D^C = B^QA, and by right-multiplieation by V, 
we get 

B'^QAV = 0 (28) 


so that 

B'^QA^V = B^QAVK^-^ =0, / = 1, 2,..., n - 1. 

Thus, for any / = 0,1,..., n — 1, we have 

B^ {A'^yu = B^ {A'^YQA^V = B^{A~^yQA^VK^-^ 
= B'^QVK^-^ = B^QAVK^-^-^ = 0 . 


(29) 


(30) 


In eonelusion, im (f/) 7 ^ {0} is eontained in the unobservable subspaee of the pair {A~^,B~^) 
and this is a eontradietion beeause {A, B) is, by assumption, reaehable, so that {A~^,B'^) is 
observable. 

6 ) This point is the dual of the previous one. 

7) Sinee P is elearly invertible we ean use the same argument employed in the proof of point 
to show that 


W := 


-Q 0 
0 


B~^ 


{-Q-YyA B] 


A'^QA - Q A'^QB 
B^QA B^QB + I 

is positive semidefinite and has rank equal to m. Therefore, there exists a full row-rank matrix 
[Co I Do] e sueh that ly = [Co | C>o]^[Co | Do]. In partieular, 

2l^Q2l-g = C^C = Co^Co 

so that there exists an orthogonal matrix U sueh that C = UCo- Let D := UDq. Therefore, 

W := [ColDoViColDo] 

= [Co I DoVU^U[Co I Z^o] = [C I D]^[C I D]. 

In eonelusion, we have 


D'^D = I + B^QB 


(31) 









11 


and 


D'^C = B'QA 


T ! 


(32) 


These two equations together with the seeond of ([^ give Q and henee, in view of point Q, 
Q{z) = C{zl — A)~^B + D is all-pass. □ 


Remark 2.1: In pointj^ of Theorem 2.1 the assumption of reaehability of (A, B) ean probably 
be eliminated for the first part of the result. More precisely, we suggest the following conjecture 
whose proof, however, seems to be non-trivial. 

Conjecture 2.1: Let A G G be given. Then, there exists P = such that 

APA~^ — P = BBA if and only if there exist matrices C G and D G such that 

Q{z) given by Q is an observable realization of an all-pass function. 

Of course, we have a dual conjecture for point [^. 

Remark 2.2: Consider an all-pass function Q{z) represented by Q. Clearly Q{z)U is still all¬ 
pass for any orthogonal matrix U. The two functions Q{z) and Q{z)U = C{zl — A)~^BU + DU 
correspond to the same dynamics so that it is natural to regard these two functions as equivalent. 
By considering the polar decomposition of D we immediately see that for any given D there is 
a unique Dq = DU such that Dq = Dj > 0. Therefore, from now on, whenever convenient, we 
can safely assume, without loss of generality, that the “D” matrix of the all-pass function Q{z) 
is symmetric and positive semidefinite. 

Remark 2.3: Consider point (or[^) of Theorem 2.1 If A is unmixed, once given A and 
B, the solution P of Q is uniquely determined and hence also the matrices C and D for which 
Q{z) = C{zl — A)~^B + D is all-pass are uniquely determined up to multiplication on the left 
side by a common orthogonal matrix. This is not the case when A is not unmixed. In this case, 
for any particular solution P of Q there exists a particular pair of matrices C and D (essentially 
different, i.e. not differing for multiplication on the left side by a common orthogonal matrix) 
for which Q{z) = C{zl — A)~^B -I- P is all-pass. Notice, however, that, once fixed A, B and 
P, the matrices C and D are always uniquely determined up to multiplication on the left side 
by a common orthogonal matrix. 

' ' '20 


Similar considerations can be made for 6). For example, let A = 
this case the set of all solutions Q of (|6) can be parametrized as Q = 


0 1/2 
1/3 


and C = J. In 


9 

-4/3 


with q 
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being a real parameter. For example, for g = 0, we get Bq = 


3 0 

0 -3/4 


and Do = 


0 1/2 

where the degree of freedom corresponding to the choice o? an arbitrary orthogonal matrix 
multiplying both Bq and Dq on the right side, has been fixed in such a way that Dq = > 0: 

since Dq is non-singular this procedure does not leave any further degree of freedom. For 


g = 1/6, we get Bi/q = 


2.85 

0.14 


0.57 

-0.71 


and Di/q = 


1.95 0.14 
0.14 0.52 


, where, again, the degree 


of freedom corresponding to the arbitrary orthogonal matrix has been fixed in such a way that 
Di/q = > 0. In conclusion, the two solutions corresponding to g = 0 and g = 1/6 lead to 

all-pass functions with different dynamical properties. 


III. LMFs AND Homogeneous Aegebraic Riccati equations 


All-pass functions can be seen as spectral factors of a spectral density function identically equal 
to the identity matrix; i.e. <I>(z) = /. This point of view turns out to be useful for classification 
of all-pass functions having a pre-assigned pole dynamics. It is a classical result in system 
and control theory lfT4l that rational spectral factorization can be cast in terms of linear matrix 
inequalities (LMI). This point of view will be used here. It will be further developed in a 
forthcoming companion paper [O devoted to stochastic modeling. In this section we shall just 
consider square spectral factors which are all-pass. 

To fix the pole dynamics we may either assign a reachable pair {A G B G or 

an observable pair {C G A G These are two “dual” structural data which will be 

fixed hereafter. Accordingly, define 


M(F) : = 


APA~^ - P APC^ 
CPA~^ CPC^ +1 


and 


N{Q) : = 


A~^QA - Q 
B^QA 


A~^QB 
B^QB + I 


and consider the two dual, constrained linear matrix inequalities (CLMI), 


M{P) > 0 
rank[M(P)] = m 

N{Q) > 0 
rank[A^(Q)] = m 


The following is an immediate corollary of Theorem 2.1 


(33) 


(34) 


r\r> A TT'X' 
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Corollary 3.1: Let P = P be a solution of ( [33] ) and let 

G 


M(F) = 

be a factorization of full rank m. Then 




(35) 


Ql{z) ■.= C{zI - A)-^G + L 


(36) 


is a (in general non minimal) realization of a square all-pass function. Dually, let Q = be a 
solution of (I34b and let 


N{Q) = 

be a factorization of full rank m. Then 


H J 


T r 


H J 


Qr{z) ■.= H{zI - A)-^B + j 

is a realization of a square all-pass function. 


(37) 


(38) 


Clearly P = 0 and Q = 0 are always solutions of the inequalities ( [33] ) and ( |34| ) and it may 
well happen that these inequalities admit no other solutions save for these trivial ones. We need 
to exclude these uninteresting circumstances. We shall henceforth assume that there is a P such 
that the matrix function with (minimal) realization 

Q{z) ■=G{zI - A)-^B + D (39) 

is all-pass. By Theorem this happens if and only equations ( [35] ) with G = B and L = D 


hold for a nonsingular P = Pq or, equivalently, if and only if ( [37] ) with H = G and J = D 
hold for a nonsingular Q = Qq. Pq and Qo turn in fact out to be such that PqQq = P In the 
next section it will be shown that each Qr{z) is a right factor of Q{z) and each Ql{z) is a left 
factor of Q{z). 


Notational convention: From now on, {A, B,G, D) such that (39) is a minimal realization of 


a square all-pass function will be the problem data; the unique solutions of (j^ and (|^ will 
be denoted by Pq and Qq, respectively and we shall reserve the symbols P and Q for generic 
solutions of ( [33] ) and ( |34| ). 


Theorem 3.1: Let (39) be a minimal realization of a square all-pass function. Then 
1) (i) For each solution P = P^ of (33), the subspace 

y = ker(P) 


(40) 
























is A^-invariant. 

(ii) The set of non-singular solutions of (33) ean be parametrized as: 


P={Pa: (41) 

where Pa := Pq is the unique solution of ([^, and Vp is the veetor spaee of 


solutions of /S.A — A = 0. If A is unmixed, then Vp = {0} and (33) admits a unique 
non-singular solution Pa = Pq, whieh is the unique solution of (|^. If A is not unmixed, 
then P eontains infinitely many solutions. 


(iii) Let Pa be a non-singular solution of (33); then to any A^-invariant subspace y there 


corresponds a solution P of ( [33] ) given by 

P:= [(/-n)p^'(/-n 


(42) 


where II is the orthogonal projector onto y. The kernel of P is y. If A is unmixed, equation 
(42), with Pa = Pq being the unique solution of (j^, parametrizes the set of all solutions 


of (^) in terms of A^-invariant subspaces. 

2) (i) For each solution Q = of (34), the subspace 


X = ker(Q). 


(43) 


is A-invariant. 


(ii) The set Q of non-singular solutions of (34) can be parametrized as: 

Q = {Qa: AePj (44) 

where Qa := (Qo ^ + A)“^, Qq is the unique solution of and Vg is the vector space 


of solutions of AAA~^ — A = 0. If A is unmixed, then Vg = {0} and (34) admits a unique 
non-singular solution Qa = Qo, which is the unique solution of (j^. If A is not unmixed, 
than Q contains infinitely many solutions. 

(iii) Let Qa be a non-singular solution of (|34|), then to any A-invariant subspace X, there 


corresponds a solution Q of ( [3^ given by 

0:= [(/- n)0i'(/- n)] 


(45) 


where II is the orthogonal projector onto X. The kernel of Q is X. If A is unmixed, 
equation ( |45| ), with Qa = Qo being the unique solution of (j^, parametrizes the set of all 
solutions of ( [34| ) in terms of A-invariant subspaces. 

Proof: We prove only point 2), as the proof of point 1) is dual. 

(i) It is clear that (34) is equivalent to existence of two matrices H G and J G such 
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that [H I J] has full row-rank and N{Q) = [H \ J]'[H \ J]. Therefore, if Q is a solution of 


(34) then A~’^QA — Q = H, so that, in view of (31 Lemma 3.1], X := ker((5) is A-invariant. 


(ii) Clearly the solution Qo of (|^ is a non-singular solution of (34) and the corresponding 


matrices H and J, introduced before, coincide with C and D of ( [39| ). Then in view of Theorem 
2 .1[ point 1^, we have 


AQ^^A'^ -Qo^ = BB'. 


-1 


(46) 


Let now Qq be another non-singular solution of (34) and Co and Dq be such that N{Qo) = 
[Co I Do]^[Co I L>o]- Equivalently, Qo is a non-singular solution of (|^ corresponding to the 
quadruple {A, B, Co, Do). Using again Theorem 2.1 point W, we have that Qq^ is a solution 


of (|^ corresponding to the same quadruple, so that, in particular, AQq^AA — = BB^. 


Comparing the latter with (46), we see that = Qq^ + A where A is a solution of the 
homogeneous Lyapunov equation AAA~^ — A = 0. If A is unmixed, this equation has a unique 
solution A = 0 so that Qo = Qo- 

Assume now that A is not unmixed. Then equation AAA~^ — A = 0 has a non-trivial vector 
space Vg of solutions and the previous argument shows that any non-singular solution Q^ of 


(34) has the form (C^o ^ It remains to show that all the elements of Q are solutions 


of (34). Observe that A[Qq^ + A]A^ — [(^o^ -I- A] = BB^ for any A G Vg. Since {A,B) is 
reachable, any Q/\ := Q^^ -f A is invertible and, in view of Theorem 2.1 point [^, there exist 
two matrices Ca and Da such that Ca{zI — A)~^B + Da is a minimal realization of a rational 
all-pass function and therefore Pa = Qa^ is the solution of (|^ corresponding to the quadruple 
{A, B,Ca, Da). This is equivalent to Qa ■= Pa^ = (Qo^ + being the solution of ([^ 


for the same quadruple so that Qa is a solution of (34) which therefore has infinitely many 
solutions. 

(iii) Let X be an A-invariant subspace. Consider an orthogonal change of basis induced by the 
matrix T = [V± | U], where the columns of V form a basis for X and the columns of V± form 
a basis for X-*“. In this basis we have 

0 


and 


T X = im 


A := T-MT = T'AT = 


Ai 

A21 


0 

A 2 


(47) 


(48) 


Partition B := T ^B = B conformably as B = 


Bi 

B 2 


. Let (5 a be a non-singular solution of 


r\r> \ TT'X' 
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(34) and let Ca and Da be sueh that N{Q/^) = [Ca \ Da]^[Ca | Da]. Equivalently, Qa is the 
non-singular solution of ([^ eorresponding to an all-pass funetion deseribed by the quadruple 
{A, B,Ca, Da)- Henee, in the new basis Qa ■= T^QaT is a non-singular solution of ^ 
eorresponding to the quadruple {A,B,Ca,Da), with Ca ■= CaT. In view of Theorem |2.l[ 
point 1^, (5 a^ is a non-singular solution of (|^ eorresponding to the same quadruple. Partition 
sueh a (5 a^ eonformably with A as 


Qa = 


11 


12 


Pl2 P 22 


(49) 


and note that it must in partieular satisfy the first equation of (|^ so that the bloek of index (1,1) 
must satisfy the redueed Lyapunov equation 


^i.Pii^i —Pii + BiPj- (50) 

Sinee the pair (^4, B) is reaehable, the pair (Ai, Bi) is reaehable as well, so that from Theorem 


2.1 point 1^, it follows that Du is invertible and there exist Ci and Di sueh that Du is the 


Q: = 


(51) 


unique solution of (|^ eorresponding to a redueed quadruple (Ai, Di, Di, Di) and henee, D{^^ 
is the unique solution of ([^ eorresponding to the same quadruple. It is now a matter of direet 
eomputation to eheek that 

' r/ 0 

0 0 

is a solution of ([^ eorresponding to the quadruple (A, D, [Ci \ 0], Di). Therefore, Q := TQT 
is a solution of ^ eorresponding to the quadruple {A, B, [Di | 0]T^, Di) and henee, it is also a 
solution of (l34l). The faet that ker[(5] = X is direet eonsequenee of ( [sT] ). We need to show that 

. I ,00 

(451 is a eoordinate-free representation of Q. By observing that T IfT = 


'T 


0 I 


and 


(j - n)g^i(j - n) = (j - n)TT' qi^tt ' (/ - n) 

Dll Di2 


= (/ - n)T 

it is a straightforward eomputation to show that 

[(/-n)Qi‘(/-n)]+ = T 


Di 2 D22 


T'(/-n), 


3-1 

11 

0 


= Q. 


(52) 


The last thing that remains to be proven is the faet that when A is unmixed, all solutions of ([34|) 


are parametrized in terms of A-invariant subspaees by (45), with Qa = Qo- We have already 


shown that in this ease (34) has a unique non-singular solution whieh eoineides with the unique 


r\r> A TT'X' 
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solution Qq of Q. The representation of the other (singular) solutions can be obtained by a 
procedure similar to the one introduced above. Indeed, assume that Q is a singular solution of 
(34) and let H, and J be such that N{Q) = [H \ J]'^[H \ J]. As already proved, ker[(5] is 


A-invariant so that we can perform a change of coordinates such that in the new basis Q has 


the structure of the right-hand side of (51), with Pn being a non-singular matrix, A has the 
structure of the right-hand side of (48), and B = [Bj \ and H = [Hi \ H 2 ] are partitioned 


conformably. It is now a matter of direct computation to check that P^^ is a solution of ([^ 
corresponding to the quadruple (Ai,Pi,Pi, J) so that Pn is a solution of (|^ corresponding to 
the same quadruple. Hence, Pn satisfies the Lyapunov equation AiPuA^ — Pn = PiP^- But 
since A is unmixed, Ai is also unmixed so that Pn is uniquely determined by the Lyapunov 
equation. As a consequence, there is a unique Q with the given kernel which necessarily coincides 
with the one given by the right-hand side of (45) with = Qo and X = ker[(5]. □ 


Remark 3.1: Let Pa and Qa denote the set of solutions of (|3^ and ([34|) described by (|42|) 


and ( |45| ) for a specific A. While when A is unmixed (and hence Pp = = {0} so that we 

necessarily have A = 0) the families Vq and Qo constitute the entire set of solutions of the 
LMI’s ([33]) and (|34|), it is not clear if this also holds for the case of a mixed A even if one 


takes the union of the sets Pa with respect to A G Pp or the union of the sets Qa with respect 
to A G Pg. The theorem provides a bijective correspondence between the family Qo of the 
solutions of ( |34| ) and the family of A-invariant subspaces. When A is not unmixed, ( |M| ), besides 
Qo, has infinitely many other families of solutions each of which being likewise parametrized 
by A-invariant subspaces. Each of these families corresponds to a non-singular solution Qa G Q 
of ( |34| ) where Q is the set of non-singular solutions parametrized by (44). The family Qo 
corresponding to A = 0 will play an important role in the following. 


Similar considerations can be made for the dual family Pq of solutions of (33) which, in case 
of unmixed A constitutes the set of all solutions of ([3^ and in case of a mixed A is just one 


of infinitely many families of solutions of (33). 

Remark 3.2: There is an obvious bijective correspondence between the set of A-invariant 
subspaces and that of A^-invariant subspaces. Indeed, X is A-invariant if and only if V := X“*“ 
is A^-invariant. This correspondence induces a bijective correspondence between the sets Pq 
and Qo. In fact, to any solution Q = [(/ — n)(5o ^(/ — H)]^ G Qo there corresponds a solution 
P = [nPQ-^n]^ G Pq. To see this, just note that the orthogonal projector Ily onto ^ := X-*- is 
equal to (/ — IT), with 11 being the orthogonal projector onto X. In this case we shall call P and 


Q complementary solutions of the LMTs (33) and (34). Indeed for complementary solutions we 


r\r> A TT'X' 
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have 


rank P + rank Q = n. 


Of course, when A is not unmixed, a similar correspondence holds for any pair of families Va 


and Qa' of solutions of ( |^ and ( |M| ) respectively, where Va is the family corresponding to a 
certain Pa G P and Qa' is the family corresponding to Qa' ■= Pa^ ^ Q (with A' ;= Pa — Pq). 


A. The case of A non-singular: Riccati equations 


In case of a non-singular A matrix, equations (33) and (34) reduce, respectively, to the 
following homogeneous algebraic Riccati equations (ARE) 


P = APA ' - APC ' (/ + CPC 




CPA 


T 


and 


Q = A^QA - A^QB{I + P^gP)'^P^QA 

The equivalence of the two representations is stated in the following proposition. 


(53) 


(54) 


Proposition 3.1: Let (39) be a minimal realization of a rational discrete-time all-pass function 
and assume that A is non-singular. Then P = P^ is a solution of (33) if and only if it is a 


solution of (53) and Q = Q is a solution of (34) if and only if it is a solution of (54). 


Proof: We prove only the equivalence of ( [34] ) and ( [54| ) as the other equivalence is dual. Let 

and J G such that N{Q) = [H \ 


g be a solution of (34). Then there exist H E 


J]'^[H I J]. In view of Theorem 2.1 pointQ, H{zl — A)~^B -f J is all-pass. After eliminating 
the non-observable part of this realization we obtain a minimal realization say C{zl — A)~^B + J 
of the same all-pass function where the A matrix clearly remains non-singular. This, in particular 
implies that J is also non-singular so that / -f B^QB = J is strictly positive definite and 
hence invertible. Then, rank[A(g)] = m implies that the Schur complement of / -f B^QB in 


N{Q) vanishes which is equivalent to Q being a solution of (54). 

Conversely, let Q = be an arbitrary solution of (54). To show that Q satisfies the LMI (34) 


it is enough to show that / -I- B^QB is positive semi-definite and, hence, positive defnite. In 
fact, in this case we can use, in the opposite direction, the previous argument based on the Schur 
complement. 


The Riccati equation (54) can be written as 


QA-^ = A^g - A~^QB{i + p^gp)-^p^g 


(55) 
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from which it is easy to see that ker(Q) is A ^-invariant and henee A-invariant. Seleet a basis 

Bi 


where A has the form shown in the right-hand side of (481, -B = 


m 




w 


is partitioned eonformably 
we Pii is non singular so that 


and Q has the same strueture of the right-hand side 

Qu ■= is also non-singular. Then substituting Q = diagjQn, 0} into (54) it is immediate 

.-1 

11 

Ff/ = AJPA^A, - AJPA^B,{I + BjPA^B,)-^BjPA^Ai (56) 


to see that P,A satisfies 


so that, by using the Sherman-Morrison-Woodbury formula, we get AiPnAl = Pn + FiF^. 
Sinee {A, B) is reaehable, {Ai, Bi) is also reaehable and Theorem 2.T point |^, implies that 
I + BjQiiBi is positive semidefinite. Observing that I + B^QB = I + BjQuBi eoneludes 
the proof. □ 


Notiee that the ARE’s (53) and (54) do not impose explieitly any positivity eondition: the 
previous result shows that these eonditions are automatieally met when A is non-singular. On 
the eontrary, when A is singular, it seems that one needs to impose explieitly the positivity 


eondition in (33) and (34): this may be merely due to a teehnieal diffieulty and we eonjeeture 


that the LMI (33) has the same solution set of the equation rank[M(F)] = m and dually, the 


LMI (34) has the same solution set of equation rank[A^(Q)] = m. 

As a direet eonsequenee of Theorem 3.1 and Proposition |3.1[ we have the following eorollary. 


Corollary 3.2: Let (39) be a minimal realization of a rational bi-proper diserete-time all-pass 
funetion. Then 

1) The unique solution Fq = Fq^ of (|^ is also a non-singular solution of the homogeneous 


Rieeati equation ( [5^ . This solution generates the family Fo of symmetrie solutions of i 
as deseribed by equation ( [42] ), where Fa = Fq and where If is the orthogonal projeetor 
onto an A^-invariant subspaee V. The elements F = F^ of this family are in a one-to-one 
eorrespondenee with the set of -invariant subspaees. 


If A is unmixed then Fq is the only non-singular solution of (53) and Fq is the set of all 


solutions of (53). 


2) The unique solution Qq = Qq of ([^ is also a non-singular solution of the homogeneous 
Rieeati equation ([54]). This solution generates the family Qq of symmetrie solutions of 


( |54] ) as deseribed by equation ( [45] ), where Qa = Qo and where 11 is the orthogonal 
projeetor onto an A-invariant subspaee X. The elements of this family are in a one-to-one 
eorrespondenee with the set of A-invariant subspaees X. 
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If A is unmixed then Qq is the only non-singular solution of (54) and Qo is the set of all 


solutions of (54). 


IV. Factorization of all-pass functions 

In this section we discuss a remarkable relation between solutions of the constrained LMFs 
(or ARE) and all pass divisors. The background facts are established in the following lemma. 


Lemma 4.1: Let (39) be a minimal realization of a rational discrete-time all-pass function 
and let Qo be the unique solution of (|^. Let P e Vq and let Q G Qo be the complementary 


solution of ( [34] ) associated to P in the sense described in Remark [3^ Let X := kerQ be the 
A-invariant subspace corresponding to Q and ^ := kerP = X-*- be the -invariant subspace 
corresponding to P. Then, one can select a basis such that, X, V, A, P, Q, Q, P and Qo have the 
following structure 


X = im 


, y = im 



Afp 

0 


Bi 

A = 



, B = 


A 21 

Tiz 


B 2 


, C* = [Qi I Qs], 


Q = 


Pn 0 
0 0 


p = 


0 0 


0 Q 


-1 

22 


) Qo — 


p-1 

11 


0 

Q 22 


(57) 


(58) 


(59) 


Proof: Perform a preliminary change of basis as in equation ( |47| ) of the proof of Theorem 
3.1] (but now use a slightly different notation) so that X, ^ are given by (57), and A has the 

Ar 0 

block-triangular structure A = 


. In this basis, partition Qo and Pq = Qo ^ as Qo = 


and Po = 


-Pii -P 12 
P 12 P 22 


A21 Ai 

We have already proved that in this basis Q = 


Qii Qi2 

Q72 Q 22 

Considering (42), where we set Pa = Qq^, and 11 is the orthogonal projector onto y, we see 

0 0 


Pn 0 
0 0 


that in the same basis P = 


0 Q-' 


22 


. Partition B and C conformably as P = 


C = [Cl \ € 2 ]- From the first of (|^ it follows that 

AJ Q 22 A 1 — Q 22 = cj C 2 - 


Bi 

B 2 


and 


(60) 


Since (A, (7) is observable, (A 2 ,Q 2 ) is observable as well so that from (60) it follows that 


Q 22 is non-singular, [|3l Lemma 3.1]. Since Q and Q 22 are non-singular the Schur complement 
Qii — Q 12 Q 22 Q 12 is also non-singular and Pn = (Qn — Qi 2Q72^Q72)~^- Perform now a further 
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change of basis induced by T = 


-inT 
12 


Q22 Q 


Although the two subspaces (57) are no 


longer orthogonal they are still in direct sum, matrix A in (581 is modified just by changing A 21 
into A 21 := A 21 + Q 2 iQj 2 ^i ~ ^ 2 Q 22 Qi 2 ’ ill we have B 2 := B 2 + Q 22 QJ 2 B 1 , and 
Cl := Cl - C2Q22QJ2. □ 

Theorem 4.1: Let ( [39] ) be a minimal realization of a rational discrete-time all-pass function. 
Let Vq be the family of solutions of ( |3^ associated to the (unique) solution Pq of (|^ and Qo 
be the family of solutions of ( [34] ) associated to the (unique) solution Qq or ([^, as described in 
Remark l3.llH 

1) For each P G Vq, G and L be such that [(7^ | has full row-rank and 


M(P) = 


'TiTr^T 


[G' 


'Ti 


Then 


Ql{z) -.= C{zI - A)-^G + L 


(61) 


(62) 


is a (non-minimal) realization of a left all-pass divisor of Q{z). The McMillan degree n; 
of Qiiz) is equal to the rank of P. 


Conversely, any left all-pass divisor of Q{z) is given by (62), where [G^ | L~^] has full 


row-rank and satisfies ( [M] ) for a suitable P e Vo- 
2) For each Q G Qo, let H and J be such that [H \ J] has full row-rank and 


Then 


N{Q) = [H I JY[H I J], 


Qr{z) -.= H{zI - A)-^B + j 


(63) 


(64) 


is a (non-minimal) realization of a right all-pass divisor of Q{z). The McMillan degree 
of Qr{z) is equal to the rank of Q. 


Conversely, any right all-pass divisor of Q{z) is given by (64), where [77 | J] has full 
row-rank and satisfies ( [63l ) for a suitable Q & Qo- 

Proof: We prove only point 1) as point 2) is dual. We first observe that Ql{z) is all-pass; 
in fact, P is a solution of ([^ associated to the quadruple (A, G, G, L) so that in view of point 
1^ of Theorem 2^ Ql{z) is all-pass. 


'As already observed, under the additional assumption that A is unmixed, Vo is the family of all symmetric solutions of 


and Qo is the family of all symmetric solutions of 134 


r\r> \ TT'X' 
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Let P e Vo and Q G Qo he complementary solutions of ( 33 ) and ( 34 ), respectively, as 


described in Remark 3.2 Select a basis as in Lemma 4.1 so that X,''^,A,B,C,Q,P and Q 


have the structure described in ( 57 ), ( 58 ) and ( 59 ). In the chosen basis, compute M{P) to obtain 


M(P) = 


0 0 

1-1 /iT ^-1 


0 


0 - Q2i 

0 C2Q22^J C 2 Q 22 CJ +/ 


( 65 ) 


0 

Gi 


so that G must have the block structure, G = 
realization 

Ql{z)=G2{zI-Ai)-^Gi + L 


and QLiz) defined in ( 62 ) has the following 


( 66 ) 


(it could be shown that this realization is minimal but we will find that this result comes as a 
byproduct at the end of the proof). Observe now that [A, G) is observable so that {Ai, G2) is ob¬ 
servable as well. Now since Q22 is a solution of (|^ associated with the quadruple {Ai, Gi,G2,L), 
then Q22 must be a solution of (|^ associated with the same quadruple. In particular, from 
the second of equations (|^ we get AJQ22G1 = GjL which may be rewritten as [AJQ22 \ 


-G 2 


Gi 


= 0 . In this factorization, the matrix [AjQ22 \ —G2] has full row-rank; in fact. 


(Ai,G^ is observable so that [Ai \ G2] has full row-rank; hence [Ai \ —G<^ has full row-rank 
lermore since Q22 is non-singular also [A^JQ22 \ —C'2] has full row-rank. The right 

has full column-rank; in fact, we have already observed that its transpose has 

full row-ranl 

dimension m and hence the m linearly independent columns of the matrix 


as well; furt 

. Gi 
matrix 


c. In conclusion, [AjQ22 \ —C'2] G has rank n so that its kernel has 

Gi 


are a basis for 


)mi[A,JQ22 I —C'2]. 

Now use the fact that Qo is a solution of equations (|^ associated with the quadruple {A, B, G, D). 
From the lower block of the second of these equations, we get 


[A^Q22 I —C'2] 


B 2 

D 


= 0 . 


( 67 ) 


Similarly, from the left-lower block of the first of the same equations, we get 

^21 
Cl 


[AjQ 


22 


-Co 


= 0 . 


( 68 ) 


r\r> \ TT'X' 
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Hence, there exist matrices and Cr such that 


B2 


Gi 

Dr, 

A21 


Gi 

Gr. 

(69) 

D 


L 


G^ 


L 




It is now a matter of direct computation to see that 

Ar 


Q{z) = [LCr I Cs] zl - 


+LDr 


0 


GlCr A; 


-1 


GiDr 


Gr 


= [G2{zI - Ai)-^Gi + L][Gr{zI - Ar)-^B^ + D^] 
= Ql{z)Qr{z) 


(70) 


where we have introduced the rational function Qr{z) := Gr{zl — Ar)~^Bi+Dr. Note that, since 
Q{z) and Ql{z) are all-pass, Qr{z) is necessarily all pass. To show that Ql{z) is a left divisor 
of Q{z) it remains only to observe that Q{z) has a minimal realization of dimension n and that 
n = ni + Ur where ni is the dimension of Ai and Ur is the dimension of Ar. As a byproduct. 


(66) is a minimal realization of Ql{z) and Gr{zl — Ar)~^Bi + Dr is a minimal realization of 
Qr{z). Finally, by construction, the McMillan degree ni of Ql{z) equals the dimension of Q 22 
or, equivalently, the rank of P. 

Conversely, let Q{z) = Ql{z)Qr{z) with Ql{z) := Gi{zl - AiY'^Bi + Du and Qr{z) : = 
Gr{zl — Ar)~^Br + Dr, being minimal realizations of all-pass functions and assume that the 
McMillan degree of Q{z) equals the sum of the McMillan degrees of Ql{z) and Qr{z). Then, 
up to a change of basis which does not affect the result that we need to establish, we have that 


A = 


Ar 0 

BiGr Ai 

G = [DiGr I Gi]; 


B = 


Br 

BiDr 


D — DiDr- 


( 71 ) 


(72) 


Hence, without loss of generality, we assume that the matrices A,B,G,D of (39) have the 


expressions (71) and (72). Since Ql{z) and Qr{z) are all-pass functions, there exist an invert¬ 
ible matrix Pi solving equations ([^ associated with the quadruple {Ai, BuGi, Di) and and an 
invertible matrix Pr solving equations (|^ associated with the quadruple {Ar, Br,Gr, Dr). By 
exploiting ( |7T] ) and ( fTT] ), it is straightforward to check that, in the selected basis, diag{Pr, Pi) 
is the (unique) solution of ([^ associated with the quadruple {A, B,G, D). Hence, we have 
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Pn = 


Pr 0 

0 Pi 


Let y = im 


be an -invariant subspaee so that 


F:= [{I-U)P,-\I-U)]^ = 


0 0 
0 Pi 


e Vo. 


(73) 


By direet computation, it is also straightforward to check that 

M{P) = I I L^] 

o"' 


(74) 


where G : = 


B, 


and L ;= Di. Now define, as in (62), the left factor Ql{,z) associated with 


P, G and L given by (73) and ( [74| ). By eliminating non-reachable part of this Ql{z), we see 
that Ql{z) = Ql{z). □ 


Theorem |4. 1 1 provides a one to one correspondence between the family Vo of solutions of (33) 
and left all-pass factors of Q{z) defined up to multiplication from the right side by a constant 
orthogonal matrix U. Similarly, Theorem |4.1|also provides a one to one correspondence between 


the family Qo of solutions of (34) and right factors of Q{z) defined up to multiplication from 
the left side by a constant orthogonal matrix U. On the other hand, a left factor Ql{z) of Q{z) 
is associated with a right factor Qr{z) by the factorization relation Q{z) = Ql{z)Qr{z). Given 
a factoriaztion of this type, it is natural to ask what is the relation between the solution P E Vo 
associated with Ql{z) and the solution Q G Qo associated with the corresponding Qr{z). The 
following result addresses this question and shows that P and Q are related by the same bijective 


correspondence introduced in Remark 3.2 


Proposition 4.1: Let (39) be a minimal realization of a rational discrete-time all-pass function 


and let Q{z) = Ql{z)Qr{z) be a minimal factorization of Q{z). The matrices P & Vo and 


Q E Qo associated with Ql(z) and Qr(z), respectively, by Theorem 4.1 satisfy the relation 
ker[P] = (ker[Q])-*- and are therefore a complementary pair. 


Proof: As in the proof of theorem 4. L let P G Pq and let Q G Qo be the corresponding 


solution of (34) as described in Remark 3.2 i.e. the only element of Qo such that ker[P] = 


(ker[Q])-‘-. We select a basis as in Lemma 4.1 so that X, V, A, P, G, Q, P and Qo have the 


structure described in (57), (58) and (59). Consider a left factor Ql{z) associated with P: as 
we have seen in the proof of Theorem 4.1 the corresponding right factor Qr{z) (that satisfies 
equation (70)) has a minimal realization of the form Qr{z) = Gr{zl — Ar)~^Bi -f Dr. Let 


Pr be the solution of (|^ associated with the quadruple {Ar, Bi,Gr, Dr). By taking (58) and 


(69) into account, we easily see by a direct computation that diag(Pr,Q 22 ) is the solution of 


r\r> \ TT'X' 
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(|^ associated to the quadruple {A, B,C, D). Since the solution P = diag(Pii, ( 522 ^) of this 
equation is unique, we have Pr = Pu. On the other hand, we know that the right factor Qr{z) 


associated with the matrix Q is given by (64) and, by duality, has a minimal realization of 
the form Qr{z) = Hr{zl — Ar)~^Bi + J. Now we compare the all-pass functions Qr{z) and 
Qr{z) and we see that they have the same state and input matrices and that the solutions of 
the equation ([^ associated to the minimal quadruple {Ar, Bi,Cr, Dr) and of the equation ([^ 
associated to the minimal quadruple {Ar, Bi, Hr, J), coincide. Hence, Qr{z) and Qr{z) differ 
for multiplication on the left side by a constant orthogonal matrix. □ 


In the case when Q{z) is bi-proper — or, equivalently, A and D are non-singular we know 
that ([^ and ([^ reduce to ARE’s. Moreover, for any given solution P of ([^, (or, equivalently. 


of (53)) we can provide an explicit expression for the matrices G and L by solving (61). The 
following corollary connects solutions of ARE’s and all-pass factorizations. 


Corollary 4.1: Eet (39) be a minimal realization of a rational bi-proper discrete-time all-pass 
function. Let Vq be the family of solutions of (53) associated with the solution P of (|^ and Qq 
be the family of solutions of (54) associated with the solution Q or (|^, as described in Corollary 

HE 


1) Eor each P G Pq, the function 

Ql{z) ■.= C{zI - A)-^G + L, 


(75) 


with 


L := (/ + GPG'^)^D 

G := APG^L-^ 


(76) 


is a (non-minimal) realization of a left all-pass divisor of Q{z). 

Conversely, any left all-pass divisor of Q(z) is given up to multiplication from the right 
side by a constant orthogonal matrix by ( [75] ), ( |7^ . 

2) Eor each Q G Qo, the function 


with 


Qr{z) ■.= H{zI-A)-^B + J, 

I J:= (J + pTQ5)i/2 

\ H := J-'^B^QA 


ill) 


(78) 


^Similarly to the general case, under the additional assumption that A is unmixed, Vo is the family of all symmetric solutions 


of 1 53 I and Qo is the family of all symmetric solutions of i54i. 


r\r> A TT'X' 
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is a (non-minimal) realization of a right all-pass divisor of Q(z). 

Conversely, any right all-pass divisor of Q{z) is given — up to multiplieation on the left 
side by a eonstant orthogonal matrix by fTT] ), ( [7^ . 

Summary and possible generalizations 

In this paper we have provided a eompletely general eharaeterization of diserete-time all-pass 
matrix funetions in the same spirit of the eontinuous-time result of Glover’s flU Theorem 5.1 ]. 
Applieations to some elass of LMI’s, to homogeneous Rieeati equations and to the faetorization 
of all-pass funetions are diseussed. The eharaeterization is presented for square all-pass matrix 
funetions but a generalization to non-square funetions ean be pursued along the same lines. 


Appendix A 

Factorization of all-pass functions which are singular at infinity 

Lemma A.l: Let Q{z) be an m x m rational proper diserete-time all-pass funetion. Then Q{z) 
ean be written as 

Q{z) = Qo{z)Qi{z)Q2{z) ... Qk{z) (79) 

where Qo{z) is a rational diserete-time all-pass funetion sueh that Qo(oo) is non-singular and 
the Qi{zys are rational proper all-pass funetions (whose only pole is in the origin) having a 
realization of the following form 


Qi{z) = 


I'm—Pi 0 


Ui + 


0 


^Pi 


(x4-o)-^[o|4]f/, 


(80) 


0 0 

where Ui is a eonstant orthogonal matrix. 

Proof: Consider a minimal realization Q{z) = C{zl — A)~^B + D. If D is non-singular, 
Qo{z) = Q{z) and we are done. If D is singular, we resort to the Silverman algorithm as 
deseribed in [31. Assume the matrix D has qi linearly independent columns, with 0 < gi < m. 


Let Vi be an orthogonal matrix such that DVi = 


Du 0 


, with Dll £ 


bmxqi 


being full column 


rank. Let us partition BVi = 


Bii B 


12 


conformably, obtaining the following block structure. 


and let 


Qi{z) := Q{z)Vi = C{zl - A) 


Qi{z) := Qi{z) 


-1 


'91 


Bqi Bt 


02 


+ 


Dm 0 


0 z T 


( 81 ) 


(82) 


r\r> \ TT'X' 
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Clearly, Qi{z) is all-pass as it is the product of all-pass functions. Moreover, Qi{z) can be 
written as 

Ql(^) = [Qll(^) I Qi2(^)] 

where 

Qii(^) Dll + CBiiz ^ -\- CABiiz ^ -f ... 


and 

^ 12 ( 2 ^) •= CB12 CAB12Z ^ CB12Z ^ . 

so that Qi{z) has the following realization 


Qi{z)=C{zI-A) 


-1 


Bii AB 


12 


+ 


Du CB 


12 


(83) 


At this point, either 


Dll CB 12 


is right-invertible, or we may iterate the above procedure 


by introducing another orthogonal matrix V 2 , such that 


Doi CB 


02 


1/0 = 


C*21 0 


with D 21 G Qf fup column rank and q 2 > qi, we define the new all-pass function 


Q2{z) := Qi{z)V2 = C{zI - A) 


-1 


B 21 B 22 


+ 


Doi 0 


(84) 


where 


B 21 B 22 


Bii AB 


12 


Co. 


Since Q{z) is all-pass, it has full rank (as a rational matrix function) and hence, after a finite 
number of steps (say k) of the above procedure, we get a rational proper all pass function 


Qk{z)=Q{z)\[V, 


i=l 


4 0 

0 


(85) 


such that Qk{oo) is non-singular. Now we set Qo{z) := Quiz), so that 


Q{z) = Qo{z) 


2=1 


4 0 

0 zlYTL—qi 


-I -1 


( 86 ) 


Finally, by setting pi := qk+i-i^ i = 1, 2,..., /c, and Ui := 14+i-i’ ^ = 1, 2,..., /c, and observing 

that ^ 

r 1 r 1 

-f 

0 0 


4 0 


0 zJ 


‘■Pi 


{zkpi 0) [0 I Ip-] 


we obtain (79) and (80). 


(87) 

□ 
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Lemma A.2: Let Q{z) be an m x m rational proper discrete-time all-pass function factored 


as in (79). Consider a reachable realization 


Qi{z) — Ci{zl — Ai)Bi -\- Di 


( 88 ) 


of Qi{z) := Qo{z)Qi{z)Q2{z) ... Qi{z). Partition Bi and A as A = [A,i | A, 2 ] and Di = 
[Di^i I 79 , 2 ], where A,i and Di i have m — pj+i columns. Then a reachable realization of 
Qi+i{z) := Qi{z)Qi+i{z) is given by 




,( 

1 

0 

0 

v 

0 / 

CJ zl - 

A,2 Aj 

J 

A,i 0 


A+1 


+[A,i I 0]A+i- 


(89) 


Proof: The realization (89) is the result of a direct computation. The fact that this realization 
is reachable may be easily seen by using the PBH test. In fact, as a consequence of the fact that 
[Ai — XI I A,i I A, 2 ] has full row-rank for all A G C, we immediately see that also 


-AJ 0 0 / 

A ,2 A ~ A,i 0 

has full row-rank for all A G C. □ 


Appendix B 

Proof of symmetry of T 
Somehow in the same spirit of we shall show that 

U := 


satisfies 


A = U-^AU, 
B = U-^B, 
C = CU. 


(90) 


(91a) 

(91b) 

(91c) 


This means that U is a similarity transform that leaves unchanged the triple (A, B, C) of the 
system. Since (A, B, C) is, by assumption, a minimal realization, this means that U = I, or that 

T = T^. 


We start with ( |91c[ ). Solving ( |12b| ) and ( |12c| ) for B we get 


B = T-^A-^C^D 


(92) 


r\r> A TT'X' 
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and 


B = AT-^C^D-^ 


By inserting in the latter the expression of D ^ obtained by transposing (12d), we get 

B = AT-^C^D - AT-^C^CA-^B 


(93) 


(94) 


Now we take the inverse of both sides of (12a) and use the Sherman-Morrison-Woodbury formula 
thus obtaining 

T-^A'^T = A-^ + A-^B{D-CA-^B)-^CA-^ 

= + A-^BD^CA-\ (95) 


From (93) we get BD^ = AT ~^C~^ whieh, plugged into the left-hand side of (95), yields 

T-^A~^T = + T-^C^CA-\ (96) 

The latter provides an expression for T~^CA~^ whieh, plugged into the left-hand side of 


(94) gives 


B = AT-^C^D + B- AT-^A'^TB (97) 

so that T~^D = T~^A^TB, or B = T~^A~^TT~^C^D. By eomparing the latter with 


(92), we eventually get 


CT ^ ^ 


-Ty^T 


(98) 


whieh, by reealling that U := T readily implies (91c). 

We now use a dual argument to obtain (91b). Solving ( |12b ) and (12c) for C we get 

C = DB'^A-^T 


and 


C = D-^B'^T'^A 


By inserting in the latter the expression of 79 ^ obtained by transposing (12d), we get 

C = DB'^T'^A - CA-^BB^T^A 


(99) 


( 100 ) 


( 101 ) 


From (|100|) we get 79^67 = B' T' A which, plugged into the left-hand side of (|95[), yields 




T ^A~^T = A ^ + A ^BB~^T~^. The latter provides an expression for A ^BB^T^ which. 


plugged into the left-hand side of (101) gives 


C = DB^T^A + C- CT-^A'^TA 


( 102 ) 
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SO that DB^T^ = CT or C = DB^T^T ^T. By comparing the latter with (99), 

we eventually get 


B^ = B'^T'^T-^ 


(103) 


whieh, by reealling that U ■.= T readily implies (|91b[). 


We now prove (91a). We multiply equation (12a) on the left side by U ^ and on the right 
side by U. By taking into aeeount ( |91b| ) and ( |91e[ ), we get 

U-^AU = T-^A-^T^ + BD-^C. (104) 


On the other hand, by transposing the first and the last member of (95) and multiplying on the 
left side by T~~^ and on the right side by we get 

A = T-^A-^T^ +T-^A-^C^DB^A-^T^ 

= T-^A-^T^ +T-^A-^C^DD-^DB^A-^T^. 


(105) 


Moreover, by inserting in the right-hand side of the latter the expressions of A and 

DB^A~^ obtained from (92) and (99), respeetively, we get 


A = T-^A-^T^ 
= T-^A-^T^ 


t/-i u 

BD-^C. 


(106) 


Finally, by eomparing the latter with (104), we get (91a). 


□ 
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